Modeling Grain Boundaries using a Phase Field Technique 
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We propose a two dimensional frame-invariant phase field model of grain impingement and 
coarsening. One dimensional analytical solutions for a stable grain boundary in a bicrystal are 
obtained, and equilibrium energies are computed. We are able to calculate the rotation rate for a 
free grain between two grains of fixed orientation. For a particular choice of functional dependencies 
in the model the grain boundary energy takes the same analytic form as the microscopic (dislocation) 
model of Read and Shockley flll| . 
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Previously we presented a new model y and @] and 
phase-field simulations of the simultaneous processes of 
solidification, impingement, and coarsening of arbitrarily 
oriented crystals. In earlier models of this phenomenon, a 
finite number of crystalline orientations are allowed with 
respect to the fixed coordinate reference frame. Morin 
et al. [[|, and Lusk Q constructed a free energy density 
having N minima by introducing a rotational (orienta- 
tion) variable in the homogeneous free energy. Chen and 
Yang J5| and Steinbach et al. ^ assigned N order param- 
eters to the N allowed orientations. In these approaches, 
the free energy density depends on the orientation of the 
crystal measured in the fixed frame — a property which is 
not physical. 

Our previous work W (KWC) introduced a model in- 
variant under rotations of the reference frame. A sim- 
ilarly motivated approach has also been developed by 
Lusk H . We have since determined that the one dimen- 
sional solutions obtained in Q (WCK) are globally un- 
stable to a lower energy solution with completely diffuse 
grain boundaries. Herein, we present a completely gen- 
eral formulation of a new class of models which allow for 
accurate, physical modeling of grain boundary formation 
by impingement and subsequent motion. 

In KWC, a rotationally invariant three-parameter 
phase field model of solidification with subsequent im- 
pingement of grains and coarsening was proposed. The 
model parameters, <fi, 77, and 9, represent: a liquid-solid 
order parameter, a coarse-grained measure of the degree 



of crystalline order, and the crystalline orientation, re- 
spectively. For what follows we shall focus on the events 
after impingement, when the material is completely solid 
(4> = 1), since we are interested in the construction of 
a physically realistic model of grain boundaries. At the 
close of this letter we will reintroduce solidification to 
complete the model. 

We now more precisely define 77 and 9. Consider a 
fixed subregion of solid material. At the atomic scale we 
can define a discrete variable which represents the ori- 
entation of an atomic bond (lattice vector) with respect 
to some fixed laboratory frame. For this subregion, we 
define 77 and 9 such that 

I N 

(77 cos 9, 77 sin 9) = — ^^(cos 6j, sinflj), (I) 

i=l 

where N is the number of bonds in the subregion. Then, 
the variable 77 is an order parameter for the degree of 
crystalline orientational order, where 77 = I signifies a 
completely oriented state and 77 = a state where no 
meaningful value of orientation exists in the subregion, 
and < 77 < 1 always holds. The variable 9 is an indica- 
tor of the mean orientation of the crystalline subregion. 
Note, that when symmetrically equivalent crystalline di- 
rections exist, 9 may not be uniquely defined. However, 
the point group symmetry operations can be used to map 
each of the 9i into a subdomain of (0, 2ir) where they are 
unique. 

Given the order parameters 77 and 9 we need to con- 
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struct a free energy which is invariant under rotations of 
the reference frame. We start from the expression 



F 



f(ri) + s/i(ri)G(\Ve\ 



ylV^I 2 



dV, 
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where v and s are constants, f{rf) is the homogeneous 
free energy and G(|V0|) is specified below. The func- 
tion u.(rj) should have the property fi(0) — and should 
increase monotonically with r\. The constants s and v 
are parameters which set the strength of the penalties 
for gradients in misorientation and degree of crystalline 
orientation, respectively. The total free energy is frame 
invariant since it includes the variable 9 only in the form 
V0. For this work we assume that only ordered phase is 
stable, so f(rj) has a single- well at rj = 1. 

In two dimensions, the equilibrium grain boundary en- 
ergy is generally a function of two angles, for example, 
the angles of a boundary with respect to the grains on 
either side. We will derive the equilibrium grain bound- 
ary energy as a function of the difference of these an- 
gles (the misorientation), A0, in one dimension. Effects 
of anisotropy (the dependence of grain boundary energy 
with respect to boundary orientation; i.e., the direction 
of V?7 with respect to mean crystalline orientation 9) are 
not currently included in our model. 

Because the energy form does not include 9 explic- 
itly, no preferred values of 9 are selected by the form 
of our equations. Therefore, values of 9 are determined 
by Dirichlet boundary conditions, if they exist, and by 
gradient flow if the value of 9 is not fixed at any spa- 
tial position. For a bicrystal (a single grain boundary) 
with Dirichlet boundary conditions, the spatial extent of 
the grain boundary is determined by G(|V0|) and the 
coupling to rj. The customary choice of square gradient, 
G(|V0|) = V0 ■ V0, cannot be physically correct, since 
it leads to a diffusion equation for 9 and therefore the 
region of misorientation spreads indefinitely. In WCK an 
attempt to stabilize a finite grain boundary was made 
by including many minima in G at different values of 
|V0|. With this form of G, a one dimensional solution 
to the equations of motion can be found which has a fi- 
nite grain boundary thickness. We demonstrate below 
that this solution has a larger energy than a completely 
diffuse interface, and is thus unstable. 

It is useful to determine the energetic contribution of 
terms with the form |V0| a ; where a is a positive num- 
ber. If the gradient is completely spread out over a one 
dimensional domain of length L and n is uniform, then 
V0 = A9/L, where A9 is the change in angle across 
the domain (grain boundary misorientation). If rj = fj 
where fj gives a minimum of f{rj) + sfj,(rf) |V0| a , and L 
is taken to be large, the total excess energy (associated 
with A9 ^ 0) is given by 



where Af(rj) = f(jf) — /(l). In particular, for a > 1 

lim AF QiL = 0. (3) 

L — >oo 

Thus, there is no form for a > 1 leading order behavior in 
G which has a globally stable, non-diffuse solution, since 
a completely diffuse boundary has an energy which can 
always be made smaller than any other configuration. 

However, the form of AF Ui l suggests that the choice of 
a = 1 may admit solutions for equilibrium grain bound- 
aries of finite spatial extent. If a < I then the derived 
evolution equation is completely singular in the region 
where V0 = 0. Thus, a = 1 (G(|V0|) = |V0| to leading 
order) is the only possible choice. 

Therefore, the dimensionless energy takes the following 
form: 
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Here we choose for f{rf) and fi(rj) the forms 



(4) 



(5) 



We examine an alternative choice for fi(n) at the close of 
this letter. 

We consider rj and 9 as a polar coordinate system on 
the unit disk D = {?; < 1} as Eqn.[l] suggests, and in- 
troduce L 2 -norm to the functional space |13|. Then, the 
equations of motion are 



= sV 



V0 



|V0| 



(6) 



(7) 



where t v and Tg are (possibly anisotropic) inverse mobil- 
ities [§. 

The form of Eqn. |?] requires careful consideration since 
it contains a singular diffusivity where V0 = 0. However, 
this singularity can be treated rigorously and controlled, 
as the work of Giga and Giga || and Kobayashi and Giga 
ifiof has indicated. A complete discussion of the dynam- 
ics of this model is reserved for a longer paper, here we 
present equilibrium one-dimensional solutions. 

The equilibrium solutions r)(x) and 9(x) must satisfy 
the following equations: 



= v 2 rjxx + 1 - V - 2sn\9 x \ 



= s 



\8x\ 



(8) 
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AF a , L = 



[Af{f 1 ) + sn{n)\V9\ a ]dx = 0{L 1 ~ a ) 1 



where the subscript x indicates differentiation with re- 
spect to x. Dirichlet boundary conditions 0(±oo) = 0± 
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are applied as well as the condition < r)(x) < 1. 
Here, let us consider the equilibrium solution which cor- 
responds to a bicrystal. Without loss of generality, the 
center of the grain boundary is located at x — 0. If f](x) 
has only one minimum, it can be shown that the solu- 
tion for 9(x) is a step function at the point where -q takes 
minimum px[ |. So we can take 9 as 



— oo < x < 0; 
< x < co. 



(10) 



Note that \9 X \ = A9 5{x), where A9 = \6+-9-\ and 8(x) 
is the Dirac delta function. Eqn. || gives 



r)(x) = 1 - (1 - i] )e » 



(11) 



where 770 — ?7(0), and should be determined by integra- 
tion of Eqn. through the discontinuity in 9: 



v 2 [VxT x Zo- = 2s m A9. 
Finally, Eqns. [ll] and [l^ determine % as 
1 „ sA6 



1 + 
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where is a scaled misorientation. 

We now calculate the total free energy. This quan- 
tity can be separated into two terms: the bulk contribu- 
tion Ebuik, and the contribution from the singular core 

at X — 0, E core . 



Ebuik = 



dx = 



e 



l + e 



and 



E C ore = S7]qA9 = V<d 



i + e 



(14) 
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With these two energies we can find the surface excess 
energy 
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Note that the energy is localized in the ^-neighborhood 
of the grain boundary, and the precise contribution to 
the energy from the region |x| > I is Ef, u ike~~ ■ 
When 



s > v, 



(17) 



then, for almost all A9, we have that 3> 1, since A9 
is expected to have the values of order one. Therefore, 
except for small values of A9, we have 



E, 



bulk 



V\ E r 



vQ- 1 < E, 



balk 



(18) 



and 



Also, 



7 — Ebulk — v. 



m ~ e- 1 < 1, 



(19) 
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which means the degree of orientational order almost van- 
ishes at the grain boundary. The typical profiles of rj and 
9 are shown in Fig [l]. Eqn. |l9| implies that the grain 
boundary energy is independent of misorientation pro- 
vided that A9 is much larger than the resolution level 
vjs which is very small in this particular case. There- 
fore, we can regard the case s 3> v to be independent of 
misorientation. 

So far, we have concentrated on the one grain bound- 
ary solution under Dirichlet boundary conditions. If the 
initial conditions had more than one boundary, the un- 
constrained grain(s) must rotate to lower the energy by 
adjusting the two boundary misorientations, until the 
equilibrium solution obtained above is attained. Since we 
are ultimately interested in the dynamics of this problem, 
it is useful to examine the case shown in Fig.^J, where a 
narrow band-shaped interior grain with thickness t lies 
between two semi-infinite grains. Like before, Dirichlet 
conditions are imposed as 0(±oo) = 9±. 

Let the orientation of the interior grain be labeled 9q, 
and assume £ 3> v, then the system is described by rj and 
9 as shown in Fig.^j. The system can eliminate the inte- 
rior grain by two independent mechanisms: 1) by moving 
the boundaries towards each other at fixed misorienta- 
tions; 2) by adjusting the orientation(s) of the interior 
grain(s) at fixed boundary spacing. Since the regions 
where rj(x) varies rapidly are spatially localized (t 3> v), 
variation of the boundary position has negligible effect. 
Therefore, we evaluate the rotation rate of the interior 
grain. If we assume, for example, 9- < 9q < 9+, the 
following estimate holds fl0| : 
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It follows from (|21j) that the angle 9q rotates to that 9± 
which is closer to 9q. If Tg is identically then the rota- 
tion rate is infinite, and a multi-grained solution would 
disappear instantly. 

As a final point, for much of this analysis we have con- 
cerned ourselves with the case where 8>1. If A9 is so 



small that <gc 1 holds, Eqn. 16 gives the asymptotics 
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sA9. 



(22) 



The energies of low-misorientation tilt grain boundaries 
have been approximated by summing the energy of dis- 
tribution of edge dislocations O]. Here we discuss a 
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method for reproducing the same energetic dependence 
on misorientation. 

It is straightforward to repeat the one dimensional 
analysis given above for arbitrary f(rj) and fj,(rj). If this 
is done for the choice 

= -2111(1-7?) -2r?, (23) 

in Eqn. ^ and leaving f(rj) unchanged, we find 77 is also 
given by Eqnjl| and t/q is determined by 

(1 - 77 ) 2 = 977 . (24) 

This guarantees < rjo < 1 for arbitrary > 0, and 
gives two asymptotic expansions of 770: 

7 l0 = 6- 1 -2®- 2 + ... /or 9»1 ! (25) 

% = 1 - 6' + . . . /or 6«1. (26) 

Thus, for 6 > 1 we obtain the same expression as given 
by Eqn. E9I In addition, for the very small values of A9 

(e « 1), 

7~-sA6»lnA(9 (27) 

holds. Eqn. ^ is the Read-Shockley energy of a low angle 
tilt grain boundary This is a very useful result, as it 
allows us to mimic the physical picture of a grain bound- 
ary as a collection of dislocations using a macroscopic 
(coarse grained) free energy for small misorientations. 

This polycrystalline model can be adopted into tradi- 
tional elements of phase field modeling of solidification; 
the full free energy will now be of the form 

F = J [/(</», 7?, t, a) + r 2 (7y, v^, e/s) (28) 

+ an{ri)\V9\ + ^\V^\ 2 ]dV. 

The field cf> is the liquid-solid order parameter, and the 
first two terms include the functional dependencies nec- 
essary to describe solidification with anisotropic kinetics 
and surface energy (see jjjj). We include the tempera- 
ture T, concentrations Cj, and any other thermodynamic 
variable which would control the state of the system. For 
anisotropic grain boundaries we need only reformulate 
j/ 2 |V7?| 2 to be functionally similar to T. With this com- 
plete model, the full process of grain formation, collision 
and coarsening can be simulated. 



To summarize, we have derived a class of models which 
can describe the solidification, impingement and coars- 
ening of grains. These models are rotationally invariant, 
and can be modified in a straightforward manner to in- 
clude a variety of physics. All of these models are ana- 
lytically tractable in one dimension, which makes study 
of their behavior particularly direct. If we choose the 
coupling carefully we are able to derive a grain boundary 
energy identical to the Read-Shockley low misorientation 
tilt grain boundary energy. Following papers will explore 
the dynamics of this model in one and two dimensions, 
and include the effects of solidification, and grain bound- 
ary anisotropy. 
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FIG. 1. Typical profiles of r/ and 9 with one grain boundary. Simulated solution and analytic solution are both drawn, which 
coincide completely within the resolution of graph. Domain size = 1.0, v — 0.01, a — 100 and A# = 
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FIG. 2. A schematic plot of the order parameter r\ and orientation 8 for a small grain with orientation 8o trapped between 
two pinned grains with orientation 8±. 
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